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Abstract 

Stripping cross sections in nitrogen have been calculated using the classical trajectory approxi- 
mation and the Born approximation of quantum mechanics for the outer shell electrons of 3.2GeV 
I~ and Cs^ ions. A large difference in cross section, up to a factor of six, calculated in quan- 
tum mechanics and classical mechanics, has been obtained. Because at such high velocities the 
Born approximation is well validated, the classical trajectory approach fails to correctly predict 
the stripping cross sections at high energies for electron orbitals with low ionization potential. 
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I. INTRODUCTION 



Ion-atom ionizine, collisions play an important role in many applications, such as heavy 
ion inertial fusion [if, coUisional and radiative processes in the Earth's upper atmosphere 
2 1, ion-beam lifetimes in accelerators [3], atomic spectroscopy and ion stopping in matter 
5|, and are also of considerable academic interest in atomic physics 0. 

To estimate the ionization and stripping rates of fast ions propagating through gas or 
plasma, the values of ion-atom ionization cross sections are necessary. In contrast to the 
electron 0] and proton ionization cross sections, where experimental data or theoretical 
calculations exist for practically any ion and atom, the knowledge of ionization cross sections 
by fast complex ions and atoms is far from complete While specific values of the cross 
sections for various pairs of projectile ions and target atoms have been measured at several 
energies [ll, l^, the scaling of cross sections with energy and target or projectile nucleus 
charge has not been experimentally mapped. 

There are several theoretical approaches to cross section calculations. These include: 
classical calculations that make use of a classical trajectory and the atomic electron velocity 
distribution functions given by quantum mechanics [this approach is frequently referred to 
as the classical trajectory Monte Carlo (CTMC) approach]; quantum mechanical calcula- 
tions based on the Born, eikonal or quasiclassical approximations, and so forth All 
approaches are computationally intensive, and the error and range of validity have to be 
assessed carefully before making any approximations or applying the results. 

Classical trajectory calculations are simpler to perform in comparison with quantum me- 
chanical calculations. Moreover, in some cases the CTMC calculations yield results close to 
the quantum mechanical calculations E3 . Es^ . The reason for similar results lies in the 
fact that the Rutherford scattering cross section is identical in both classical and quantum 
mechanical derivations |16]. Therefore, when an ionizing collision is predominantly a conse- 
quence of electron scattering at small impact parameters close to the nucleus, the quantum 
mechanical uncertainty in the scattering angle is small coinpared with the angle itself, and 



the classical calculation can yield an accurate description [17|, [ISj . But this is not always a 
case, as we demonstrate below. For fast projectile velocities and low ionization potentials, 
the difference between the classical and quantum mechanical calculations of ionization cross 
section can be as large as a factor of six for parameters to relevant to heavy ion fusion cross 
sections. 

In the present analysis, we consider at first only the stripping cross section of loosely 
bound electron orbitals of /~ and Cs~^ ions colliding with a neutral atom of nitrogen, or 
with a fully stripped nitrogen ion with Zt = 7 (for comparison). Atomic units are used 
throughout this paper with e = h = rrie = 1, which corresponds to length normalized to 
flo = h'^/irriee'^) = 0.529 ■ 10~®cm, velocity normalized to vq = e^/h = 2.19 ■ lO^cm/s, and 
energy normalized to Eq = rrieVQ = 2Ry = 27.2eV, where Ry is the Rydberg energy. The 
normalizing coefficients are retained in all equations for robust application of the formulas. 
For efficient manipulation of the formulas, it is worth noting that the normalized veloc- 
ity is v/vq = 0.2^y E[keV/ amu], where E is energy per nucleon in keV/amu. Therefore, 
25keV/amu corresponds to the atomic velocity scale. 

The typical scale for the electron orbital velocity with ionization potential /„/ is Vni = 
Vq ^y2Ini/EQ. Here, n, / is the standard notation for the main quantum number and the 
orbital angular momentum quantum number [l6|. The collision dynamics is very different 
depending on whether v is smaller or larger than Vni- 
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FIG. 1: Shown in the figure is a comparison of the ionization probabiUties [PqupiQ) in Eq.Q, and 
Pdp{q) in Eq.©] and the effective charges [ZquT^q) Nexiq) in Eq.((7|, and Zciriq) in Ea. lfTn]) ] in 
quantum and classical mechanics for 3.2GeV I~ ions colliding with a nitrogen atom. Ionization of 
only the outer electron shell is considered (here, Inip = 3eV). 



II. BEHAVIOR OF CROSS SECTIONS AT LARGE VALUES OF PROJECTILE 
VELOCITY V > Vni 

When V » Vni, the projectile interaction with the target atom occurs for a very short 
time, and the interaction time decreases as the velocity increases. For 3.2GeV I~ ions, 
envisioned for heavy ion fusion applications, the projectile velocity in atomic units is 32vq, 
while the electron orbital velocity is Vni = 0.5vq for the first (3.06eV^) ionization potential of 
/~, and Vni = 1.3fo for the first {22AeV) ionization potential of Cs~^. Therefore, we shall 
use the limit v >> Vni- 

In the limit, where v > vqZt and v » Vni, the Born approximation of quantum me- 
chanics can be used Ullil. The first inequality assures that the nitrogen atomic potential 



can be taken into account as a small perturbation (the Born approximation); the second 
inequality allows us to use the unperturbed atomic wave function. 

In both classical mechanics and in the Born approximation, the ionization cross section 
can be recast in the form l^, 3, 3 2^ - 



where Pp{q) is the probability of electron stripping from the projectile when the electron 
acquires the momentum q, and da/dq is the differential cross section for scattering with 
momentum q. 
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FIG. 2: Plots of differential cross sections for stripping of I ions by nitrogen atoms and fully 
stripped ions. 



In quantum mechanics, Pqupil) can be expressed by the square of the corresponding 
matrix element of transition from the initial state \nl > to the state of the ejected electron 
|k > with momentum k, integrated over all k. This gives 

Pgup{q) = j |<n/|e*'ink>|'rf3k. (2) 

The analytical form of Pqup{<l) for hydrogen-like electron functions is given in Ref. 19]. In 
classical mechanics, Pcip{q) is given by the integral over the electron velocity distribution 
function /(vg) defined by 

Pcip{q) = 1 e (^q ■ V, + ^ - /(Ve)dv^. (3) 

Classical mechanics prescribes the electron velocity distribution function (EVDF) for 
hydrogen-like orbitals as a microcanonical ensemble, where 



Here, C is a normalization constant defined so that / f (ve) dv^ = 1, and denotes 
the Dirac delta-function. Interestingly, the EVDF for a hydrogen-like electron orbitals is 
identical in both the quantum mechanical and classical calculations fli^, with 



/(-e)=^-^, (4) 
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where Vni is the scale of the electron orbital velocity defined by 



Vnl = Vo^/2Inl/Eo. 

In the Born approximation of quantum mechanics, da/dq is given by 3, 21 



da 
dq 



ina, 



where 



^quT 



nl 



NeT{q) = [Nir-J2\^niT{q)\' 



(5) 



(6) 



(7) 



nlT 



Here, Zqux^q) is the effective charge, subscript qu stands for quantum mechanics, Fniriq) = 
J e^"^'^ pniT{f)d^T is the form factor of the target atom's orbital nl with the electron density 
Pnirir), and N^^"''- is the total number of electrons in the target atom [Neriq oo) = 

In classical mechanics, da/dq is given by 



TVftotal 



da 2 

dq dq 



(8) 



Here, p{q) is the impact parameter for a collision resulting in the momentum transfer q. For 
fast collisions, q is mainly perpendicular to the projectile velocity, and q is determined by 
integration of the electric field of the target atom on the electron, which gives 



,(p) ^ r 

V Jo 



dUj 



dr ^ 



-dr. 



(9) 



where Uxir) is the atomic potential of the target atom. To compare the classical calculation 
with the quantum mechanical calculation, we recast Eqs.® and (jH)) into a form similar to 
Eq.®, introducing the effective charge Zcixiq) defined by 



Zcixiq) 



qv 



dp 



(10) 



where subscript cl stands for classical mechanics. Note that for the bare target ion, Ut = 
—e^Zrjr and Zciriq) = Zt- Finally, making use of the effective charge in Eq. pOj) . the 
differential cross section in classical mechanics takes on a form similar to Eq.® in quantum 
mechanics, i.e., 

da „ 2 VoimeVof ZdT{q? + K 

^0 „.2 



Ttotal 



(11) 



dq " q' 

Here, the final term accounts for ionization by the 1^^^°^ target electrons. 

Figure 1 shows a comparison of the ionization probabilities \Pqup{(l) in Eq.Q, and Pdp{q) 
in Eq.Q] and the effective charges [Zguriq) in Eq.©, and Zdriq) in Eq.ijTUl)] in quantum 
mechanics and classical mechanics for 3.2GeV I~ ions colliding with a nitrogen atom. Ion- 
ization of only the outer electron shell is considered (here, Inip = 3.06eV^, approximating as 
a hydrogen-like orbital). 
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Figure 2 shows that for stripping by neutral atoms, the main contributions arise from in- 
termediate momenta in the range q = 0.5 — 1, while for stripping by the bare target nucleus, 
small values of q make the largest contribution to the cross section, which corresponds to 
large impact parameters (due to the Coulomb long-range interaction). Because PguP > PdP 
for g << 1, but Zg^r < Zdriq), the quantum mechanical cross sections are larger than the 
classical stripping cross sections for stripping by the bare nucleus, but smaller than the clas- 
sical stripping cross sections for the atoms. Carrying out the integration in Eq. gives the 
stripping cross sections for only one electron from the outer electron shell for different ions 
with the same velocity v = 32vo colliding with a nitrogen atom. The results are shown in 
Table 1 for 3.2GeV I' ions; in Table 2 for 3.35GeV Cs+ ions; and in Table 3 for 25MeV R- . 

cr, 10~^^cm^ quantum classical 
N 0.08 0.47 

2.5 1.29 

Table 1. Cross section for stripping of 3.2GeV I~ ions colliding with a nitrogen atom 
and a fully stripped nitrogen ion ( stripping of only one electron from the outer electron 
shell is considered here with Inip = 3.06eV^ ). 

cr, 10~^^cm^ quantum classical 
N 0.045 0.10 

0.32 0.17 

Table 2. Cross section for stripping of 3.35GeV Cs^ ions (the same velocity as 3.2GeV 
I^) colliding with a nitrogen atom or a fully stripped nitrogen ion ( stripping of only one 
electron from the outer electron shell is considered here with Inip = 22AeV ). 

cr, lO^^^cm^ quantum classical 
N 0.10 1.34 

12.5 5.05 

Table 3. Cross section for stripping of 25MeV^ H~ ions (the same velocity as 3.2GeV 
I~) colliding with a nitrogen atom or a fully stripped nitrogen ion ( stripping of only one 
electron from the outer electron shell is considered here with Inip = 0.75eV ). 

Figure 3 shows the same results as in Fig. 2, but the results are obtained for 3.35GeV^ Cs~^ 
ions (ionization of only one outer electron shell is considered here with Inip = 22AeV ). 
Note that 3.3bGeV Cs^ is chosen to have the same velocity as a 3.2GeV I~ ion. 

In the limit v » Vnu the stripping cross section by a fully stripped ion can be analytically 
evaluated. The Bohr formula, derived by means of classical mechanics, neglects the electron 
atomic velocity, and gives for the cross section 

Accounting for the electron atomic velocity gives an additional factor of 5/3 jl5j]. The Bethe 
formula il9|J derived by means of the Born approximation of quantum mechanics gives 



a = cr [v.lnuZp) 



0.566 In ( — ) + 1.261 

\Vnl 



(13) 
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FIG. 3: Plots of the differential cross sections of ionization for Cs"*" and / ions by nitrogen atoms 
and fully stripped ions. 

The results of cross sections calculations using Eq. p2|l with a factor 5/3 and the result 
in Eq. p3|) coincide with the results in Tables 1, 2 and 3 of stripping cross sections by a 
fully stripped nitrogen ions calculated in classical trajectory approximation and the Born 
approximation of quantum mechanics, respectively. 

The stripping cross sections calculated in classical trajectory approximation for Cs^ and 
/~ ions by fully stripped nitrogen ions is only factor 2-3 larger than the stripping cross 
sections by neutral nitrogen atoms, which is in qualitative agreement with the observations 
in Ref. |l2j|. However, there is a large difference, up to a factor 30, in the stripping cross 
sections calculated in the Born approximation of quantum mechanics. 

It is evident that the stripping of Cs^ ions by fully stripped nitrogen ions decreases by 
a factor of 22AeV /?)eV = 7.5 compared with I~ ions, which is in agreement with the Bohr 
[Eq. ()12j) ] and Bethe [Eq.()13|)] formulas. The stripping cross sections for Cs~^ and /"ions 
by neutral nitrogen atoms differ by only a factor of 2. In classical mechanics, because the 
interaction potential is a strong function of the separation, to transfer a considerably larger 
momentum requires a rather small decrease in impact parameter. This is why, notwith- 
standing the large difference in ionization potential by a factor of 7, the difference between 
the two cross sections is only a factor of 2. Table 3 shows that the difference between the 
quantum and classical treatments increases for smaller ionization potentials (compare Table 
3 with Table 1). 

The reason for such a large difference between the quantum mechanical and classical 
mechanical stripping cross sections for /~ can be easily understood from the example of 
elastic electron scattering from the shielded Coulomb potential U{r) = exp(— r/ao)/r. The 
differential cross section for elastic scattering is shown in Fig. 4 . For the shielded Coulomb 
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FIG. 4: Plots of the differential cross sections for the shielded Coulomb potential for v = 32vq. 
potential, direct application of the Born approximation gives ^ 



qdq ° {q' + minyalY' ^ ' 

and the total cross section is a = A'ko^Vq/v'^ . The total classical cross section, obtained from 
integrating j pdp, diverges because of the contributions from large p (small q). Evidently, 
the quantum mechanical cross section departs from the Rutherford scattering formula for 
q/{meVQ) < 1, whereas the classical mechanical cross section departs from the Rutherford 
scattering formula only for q/{meVo) < 2vo/v [see Eq.(jni) and Fig. 4]. Therefore, the classical 
differential cross section differs from the quantum mechanical result by a factor of [v/{2vq)]^, 
which for v = 32f o gives a difference in small-angle differential cross section of up to a factor 
of 10^ (see Fig.4). 

Tables 4 and 5 are similar to Tables 1 and 2, but the calculations are carried out for ion 
energies 30 times smaller, in the range of lOOMeV. Table 5 shows that the predictions of the 
classical and quantum mechanical theories are similar for lOOMeV ions. However, they are 
a factor two different for /~ ions, and the cross sections are the same within 10% accuracy 
for Cs^ ions. The contribution from small q to the stripping cross section by a neutral 
nitrogen atom is smaller for Cs^ ions than for ions, thereby significantly reducing the 
stripping cross section of Cs^ ions compared with /~ ions, especially for the calculation in 
the classical trajectory approximation (see Tables 4 and 5, and Fig. 5). 

cr, 10~^^cm^ quantum classical 
N 2.47 6.8 

61 37 

Table 4. Cross section for the stripping of 105MeV /" ions {v = 5.75fo) colliding with 
a nitrogen atom and a fully stripped nitrogen ion (stripping of only one electron from the 
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FIG. 5: Plots of the differential cross sections for stripping of lOOMe VCs"*" and lOSMeV /~ ions 
{v = 7.5vo) by nitrogen atoms. 

outer electron shell is considered here with Inip = 3eV ). 

(7, 10~^^cm^ quantum classical 
N 1.36 1.4 

6.6 5.2 

Table 5. Cross section for the stripping of llOMeV Cs^ ions {v = S.TSwq) colliding with 
a nitrogen atom and a fully stripped nitrogen ion (stripping of only one electron from the 
outer electron shell is considered here with Inip — 22AeV). 

III. CALCULATION OF TOTAL CROSS SECTIONS 

The total stripping cross section is defined as 

a*<'*»' = J]m(j^, (15) 

m 

where am is the cross section for stripping m electrons in each collision. This cross section 
is convenient to use for electron production calculations. The stripping cross section for any 
degree of ionization is defined as 

a = ^am: (16) 

m 

which is a convenient expression to use to determine the ion confinement time in an accel- 
erator. In the limit v » Vni, the calculation of the total stripping cross section can be 
performed assuming that the stripping from different electron orbitals occurs independently 
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[i3,i.e., 

(y'"'"' = Y,^ni(yni, (17) 

nl 

where ani is the stripping cross section of only one electron from the electron orbital nl, 
and Nni is the number of electrons in the orbital. The structure of the electron orbitals for 
ions is shown in Table 6. 

nl 5p 5s 4d 4p 4s 3d 3p 3s 2p 2s 

N„i 6 2 10 6 2 10 6 2 6 2 

Ini 3.08 13.2 50.1 125.0 185.83 623.26 892.5 1.07e3 4.65e3 5.2e3 

ani{v = 32vo) 0.080 0.054 0.030 0.018 0.013 5.5e-3 4.2e-3 3.6e-3 8.3e-4 7.3e-4 
ani{v = 5.75vo) 2A5 1.65 0.92 0.52 0.39 0.12 0.078 0.062 5.8e-3 4.6e-3 

Table 6. The structure of electron orbitals for /~ ions and the individual cross sections 
avaluated for an orbital electron in units of 10~^^cm^. 

Here, nl denotes the atomic orbital quantum numbers, Ini is the ionization potential in 
eV, and ani denotes the individual cross section for an orbital electron in units of 10~^^cm^. 
The sum over all orbitals gives a*"*"' = 1.1 • 10~^^cm^ for 3.2GeV I~ ions. To correctly 
account for multiple ionization, the inclusion of multi-electron effects is necessary. This will 
be addressed in a future publication. However, it is clear that the stripping cross section for 
any degree of ionization by neutral atoms is limited by the geometrical cross section of the 
atom (the geometrical cross section of a nitrogen atom is much smaller than the geometrical 
cross section of a Cs"*" ion or a /~ion The nitrogen atom geometric cross section is 

ctat = 1.5 ■ 10~^^cm^|22], and therefore a < uat is expected. Preliminary estimates suggest 
that single electron stripping is expected under these conditions. 

For 105MeV /~ ions, however, the sum over all orbitals gives a*"*"' = 33 ■ 10~^^cm^, 
whereas ctjv = 1-5 ■ 10~^^cm^. This indicates that multi-electron ionization is expected. 
However, it is clear that the stripping cross section for any degree of ionization is limited 
from above by = 1.5 ■ lO^^^cm^. 

The structure of the electron orbitals for Cs~^ ions, and the individual cross sections for 
an orbital electron in units of 10~^^cm^ are illustrated in Table 7. Note that a Cs~^ ion has 
the same number of electrons on each orbital as a /~ ion. 

nl 5p 5s 4d 4p 4s 3d 3p 3s 2p 2s 

N„; 6 2 10 6 2 10 6 2 6 2 

Ini 22.4 34.0 88.3 176 242 742 1.03e3 1.2e3 5.1e3 5.7e3 

ani{v = 32vo) 0.044 0.037 0.022 0.014 0.011 4.8e-3 3.7e-3 3.2e-3 7.4e-4 6.5e-4 

a^ilv = 5.75vo) 1.35 1.12 0.66 0.41 0.32 0.098 0.065 0.052 4.7e-3 3.8e-3 

Table 7. The structure of electron orbitals for Cs^ ions and the individual cross sections 
for an orbital electron in units of 10~^^cm^. 

For 3.35GeV Cs~^ ions colliding with a nitrogen atom with velocity v = 32vo 
{25MeV/amu), the summation in Eq. (jl7p over all orbitals gives cr*°*'^' = 0.72- 10~^^cm^. This 
estimate of the cross section is consistent with Olson's result in Ref. [12], a = 2 ■ 10~^^cm^ 
for 25MeV/amu Xe^ . Note that the factor of three difference between the results presented 
in Table 7 and the results in Ref.llJl is due to the fact that the cross sections in Table 7 
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are predicted by making use of quantum mechanics, whereas results in Ref. are classical 
trajectory calculations, not applicable at such high projectile velocities. 

For llOMeV Cs~^ ions colliding with a nitrogen atom, v = 5.75vo {O.SMev/amu) and the 
summation over all orbitals in Eq.(|17p gives cr*"**^' = 21 ■ 10~^^cm^, whereas the geometrical 
cross section of a nitrogen atom is only cttv = 1-5 • lO^^^cm^ << cr*"*"'. This indicates 
that multi-electron ionization is expected, similar to ions at the same velocity. As noted 
earlier, to correctly account for multiple ionization, multi-electron calculations are necessary. 
However, it is clear that the stripping cross section a for any degree ionization is limited 
by CT/v = 1.5 ■ 10~^^cm^. This estimate of the cross section is consistent with Olson's result 
l^, £7*°*"' = 4 ■ 10"^^cm^ for 2MeV/amu Xe+. The inequahty cr*°*"' > ctn indicates the 



important effect of multi-electron events. 



IV. CONCLUSIONS 

For low ionization potential, where a small momentum transfer q contributes to stripping, 
the classical approach is not valid. For 3.26*6^^ /~ ions, the classical trajectory approach 
overestimates by a factor of six the stripping cross section by atomic nitrogen, and by a 
factor of two the stripping cross section of 3.35GeV Cs^ ions. For llQMeV Cs^ ions and 
IQbMeV I~ ions colliding with a nitrogen atom at velocity v = 5.75vq {0.8Mev/amu), 
multi-electron ionization is expected. For a correct description of multiple ionization, multi- 
electron calculations are necessary. However, it is clear that the stripping cross section for 
any degree of ionization is limited from above by the geometrical cross section of nitrogen, 
with (Tat = 1.5 ■ 10~^^cm^, and should be be similar in magnitude for /~ ions and Cs~^ ions 
at energies in the lOOMeV range. (The geometrical cross section of a nitrogen atom is much 
smaller than the geometrical cross section of a Cs^ ion or a J~ion i^]. This effect is similar 



to the hole produced by a bullet piercing a paper target, where the hole size is determined 
by the bullet cross section, not by the paper target.) 
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